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ABSTRACT 

We present Gemini GMOS longslit spectroscopy of the isolated SO galaxy NGC 3115. 
We have determined kincmatical data and Lick/IDS absorption line-strength indices 
for the major axis out to around 9 kpc and for the minor axis out to around 5 kpc 
(around 2R e ). Using stellar population models which include the effects of variable 
[a/Fe] ratios we derive metallicities, abundance ratios and ages for the stellar popula- 
tion of NGC 3115. We find that [a/Fe] remains fairly constant with increasing radius 
at around [a/Fe] = 0.17 for the major axis but increases rapidly for the minor axis to 
around [a/Fe] = 0.3. We also find that to first order this behaviour can be explained 
by a simple spheroid + disc model, where the spheroid has [a/Fe] = 0.3 and the disc 
shows close to solar abundance ratios. The disc also appears considerably younger than 
the spheroid, having an age of around 6 Gyr compared to 12 Gyr for the spheroid. We 
compare these results to those previously presented for the globular cluster system of 
NGC 3115. 

Key words: galaxies: general - galaxies: abundances - galaxies: individual: NGC 
3115 - galaxies: stellar content - globular clusters: general 



1 INTRODUCTION 

Globular clusters (GC's) are among the simplest stellar sys- 
tems. The stellar content of individual clusters form primar- 
ily at one epoch and location and so are remarkably uniform 
in terms of metallicity, age and chemical abundances. Be- 
cause GC's are believed to form preferentially during the 
major star forming and mass accumulation epochs, a care- 
ful examination of the GC population of a galaxy can shed 
light on these periods of a galaxies' development. This how- 
ever relies on the assumption that GC's act as good tracers 
of the properties of the overall stellar population formed at 
the same epoch. For a more th orough discussion o f why this 
is believed to be the case see iPuzia et "all J2005). It is the 
aim of the present work to test this hypothesis in detail for 
one well studied galaxy. 

NGC 3115 is one of the closest and most stud- 
ied SO galaxies. Its GC system has been extensively 
investigated, using both photometric llKavelaarsl Il998t 
IKundu fc Whitmord Il998l: Puzia et all 120021) and spec- 



troscopic~ lKmu^clmer^itHd" ]2^0?1Pu^iaet^Ll 120041) tech- 
niques. These studies find two GC sub populations of 
mean metallicitie s [Fe/ Hl ~ -0.37 and [Fe/H] ~ -1.36. 
iKuntschner et all {2002) find that both GC sub popula- 
tions have ages which are consistent with a single epoch 
of formation about 12 Gyr ago. The observations of 
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IKundu fc Whitmord <199Sl) and iKavelaarsI <199Sft are con- 
sistent with the metal-rich clusters being associated with a 
rapidly rotating thick disc system and the slower rotating 
metal-poor clusters being associated with the halo of NGC 
3115. The stellar populatio n of the galaxy its elf has also 
been examined, notably by iFisher et al.l lll996l) who mea- 
sured the line strengths and their gradients out to a radius 
of 40" along the major axis, iTrager et al.l (|l998) also exam- 
ined t he absorption line indices of the galaxy nucleus. lElsonl 
lll997l) found evidence for a bimodal distribution of metal- 
licities with [Fe/H] ~-0.7 and [Fe/H] ~-1.3 in the resolved 
stellar population of the halo located 8.5' east of the centre 
and 5 ' from the major axis, although IKundu fc Whitmord 
( 1998) suggested that this might be due to an instrumental 
effect in the metallicity calibration. Taken together these 
observations lead to the possibility that the metal poor sub- 
population formed 12 - 13 Gyr ago, with the metal rich 
population forming several Gyr later after the ISM has been 
enriched by a factor of ~4, perhaps triggered by AGN ac- 
tivity or a merger event. A further examination of the ages 
and metallicities of the general stellar population and the 
GC population in this galaxy using updated Simple Stellar 
Population (SSP) models can test the plausibility of these 



scenarios. 



This paper is structured as follows. Section 2 describes 
the data set and its reduction. Section 3 discusses the kine- 
matical results of this investigation, Section 4 examines 
the line indices, and their use for calculation of luminos- 
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Table 1. NGC 3115 Basic Parameters 



Parameter 



Value 



Source 



Right Ascension (J2000) 
Declination (J2000) 
Morphological Type 
Major Diameter 
Minor Diameter 
Heliocentric Radial Velocity 
Asymptotic Radial Velocity 
Central Velocity Dispersion 
Inclination of disc 



10 ,1 05 m 13.98 s 
-07°43'06."9 

scr 

7.2 arcmin 
2.5 arcmin 
663 ±4 kms" 1 
263 ±5 kms" 1 
314 ±4 kms" 1 
86° 



a NED. http://nedwww.ipac.caltech.edu/ 

6 This study. 

c Capaccjrjli^^alj (1993). 

ity weighted abundance ratios, age and metallicity. Section 
5 describes a simple two component fit to the [a/Fe] ratios 
of NGC 3115, Section 6 includes a general discussion of the 
results of the previous sections and Section 7 presents our 
conclusions. 



2 OBSERVATIONS AND DATA REDUCTION 

The observations were carried out on t he 18/19th of De - 
cember 2001 with the CMOS instrument llHook et alJl2004h 
on the Gemini North telescope (Program ID GN-2001B-Q- 
44). The B600 grism with 600 lines/mm was used with a 
longslit 1 arcsec wide by 335 arcsecs long. The data were 
binned by 4 in the spatial dimension and 2 in the spectral 
dimension producing a spectral resolution of ~4.4 A FWHM 
(110 kms -1 ) sampled at 0.9 A/pixel. The seeing throughout 
the observations was generally 0.8 arcsecs and the binned 
pixel scale was 0.3 arcsecs/pixel. The wavelength range is 
~3800-6400 A. Two sets of integrations were completed, one 
each for the major and minor axis of the galaxy, with a total 
integration time for each axis of 7200s. The centre of the ma- 
jor axis longslit was offset from the galaxy centre along the 
major axis by 120 arcsecs as can be seen in Fig. 1, to extend 
the radial coverage and reduce the galaxy light in the sky 
background region. For the minor axis the sky region was 
defined at either end of the longslit. As can be seen from Fig. 
1 the sky subtraction regions are located at large radii where 
the galaxy flux is only 5 - 6% of the sky level, hence over- 
subtraction of galaxy light should not be a major problem. 
This is discussed in more detail in Sect. 2.1. Two velocity 
standards HD97907, HD73665 and a photometric standard 
Feige 34 were also observed using the same experimental set 
up. 

The standard Gemini IRAF routines were used to carry 
out bias subtraction, flat-fielding, and cosmic ray subtrac- 
tion. In addition to these procedures it was noted that scat- 
tered light seriously compromised the measurement of line 
indices at large galactocentric radii. The scattered light had 
the effect of decreasing measured line indices by applying 
a constant offset to the spectra; this effect leads to the 
measurement of spurious age gradients at large radii. Fortu- 
nately the GMOS set-up utilised in this investigation pro- 
vided the opportunity to remove this scattered light. This 
was possible because the 2-D spectrum produced contained 

3 unexposed regions created by the bottom of the slit and 
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Figure 1. Positions of the major and minor axis longslits. Slit 
widths exaggerated by a factor of 5. Dark regions define sky- 
subtrac tion regions. Filled circles denote positions of GC's exam- 
ined bv iKuntschner et alJ (2002) with line ind ex measurements, 
unfille d circles have only kinematical data fromfKuntschnc r et alJ 
(2002) available. The ellipse indicates the position and orienta- 
tion of the galaxy and enc loses roughly half the integrated light 
iMichard fc Marchallll994^ . 



two slit spacers located at 1/3 and 2/3 of the distance along 
the spatial dimension of the image. After bias subtraction 
these regions should contain no flux, however light scattered 
within the spectrograph meant that this was not the case. 
To correct for this effect it was possible to interpolate the 
scattered light level between the 3 unexposed regions and 
to subtract this from the image. In total scattered light ac- 
counted for a ~ few % of the total incident light. This how- 
ever was sufficient to account for ~ 50% of the counts at the 
largest radii. After correction residual scattered light should 
be of a negligible level. 

The data was then wavelength calibrated with the wave- 
length calibration being accurate to ^ 0.2 A. The 2-D spec- 
trum was then extracted into a series of 1-D spectra which 
were sky subtracted and binned in the spatial dimension 
until the target S/N (measured in a region near to the H/3 
line) of 20 or 60 (see Sect. 3) was reached. After binning to 
S/N = 60 a total of 108 (major axis) and 44 (minor axis) 
spectra were produced, it is from these spectra that all of 
the kinematics and line indices were measured, though in 
later figures these have been rebinned in groups of 4 (radi- 
ally) to allow clearer presentation. The binned 1-D spectra 
were then flux calibrated using the photmetric st andard and 
a red dening correction of E(B - V) = 0.146 iSchleeel et alJ 
1998) applied using the standard IRAF routines. 
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Figure 2. Typical rest-frame spectra corrected to the Lick/IDS resolution obtained in this study. Spectrum A is from the central region, 
Spectra B and C are the largest radii spectra from the major (f» 130 arcsecs) and minor axes (si 70 arcsecs) respectively. Over-plotted 
are the index bands (light grey) and red/blue continua (dark grey) for the main Lick indices measured in this study. 



2.1 Sky Subtraction 

Clearly an accurate sky subtraction is a necessity when at- 
tempting to measure accurate kinematics and line indices 
at large galactocentric radii. If the sky spectrum is strongly 
contaminated by diffuse stellar light it can introduce large 
errors in measured indices and kinematic quantities such 
as a, I13 and h.4. To investigate the magnitude of this ef- 
fect several methods were used, the simplest of which was 
to re-reduce our data using the minor axis sky to sky sub- 
tract the major axis and vice versa. The line indices were 
remeasured and compared to those measured previously, for 
all indices considered here the changes in measured index 
were < 0.1 A except for the Mgfo line which showed quite 
considerable variation. This change in the Mgb line can be 
understood as being entirely due to changes in the strength 
of the 5200 A sky emission feature between the major and 
minor axis exposures. 

The second test carried out was to reduce the data as- 
suming ± 20% errors in the sky spectra, again any changes 
in the measured indices were < 0.1 A for all cases except for 
the very last (and largest) radial bin. We are therefore con- 
fident that errors introduced by the sky subtraction should 
be entirely negligible for all but the very largest (and most 
binned) radial bins. 



2.2 Kinematics 

The first step in analysing the spectra was to measure the 
line of sight velocity distribution (LOSVD) for each of the 
binned spectra. The publicly availa ble IDL implementation 
of the pPXF method iCappellari fc Emselleml2004) was em- 
ployed to determine the recessional velocity, velocity disper- 
sion, and the Gauss-Hermite moments /13 and /14 for the indi- 
vidual spectra. In this method the NGC 3115 spectra were fit 
in pixel space over the wavelength range 4900 - 5450 A using 
both the stellar velocity templates observed during the run 



and ad ditional stellar templates from the library of lVazdekisI 
(1999). Errors were estimated by remeasuring 100 Monte- 
Carlo simulations of the spectra with added photon noise. 
Additionally, as a check of accuracy, the redshift and veloc- 
ity dispersion were measured by Fourier cross-correlating the 
spectra against the velocity standards using the fxcor im- 
plementation in IRAF. As the two methods produce consis- 
tent results, the more comprehensive pPXF implementation 
was used for the remaining analysis. 

2.3 Line Indices 

The system of line strength indices most commonly 
used for optical studies is the LICK/ID S system 
JWorthev fc Ottavianil II99I: iTraeer et af]ll99Sl) . We con- 
volve our spectra with a wavelength dependent Gaussian 
kernel to reproduce the LICK/IDS resolution of ~9 - 11 A. 

The line strength indices were then measured for each 
spectrum and corrected for the smearing effects of the 
LOSVD following the procedure of lKuntschiierl (120041) . This 
procedure corrects the line index measurements for the ef- 
fects o f velocity dispersio n and the non-gaussian /13 and hi 
terms. iKuntschneil |.2004) finds that for changes of ±0.1 in 
hi with constant a = 250 kms -1 the LOSVD correction 
changes by ±5% with corresponding errors in the age and 
metallicity estimations of 15 - 20%. Clearly for a galaxy such 
as NGC 3115, with central velocity dispersion of around 300 
kms -1 , the effect of these corrections on the spectra from 
the central few arcsecs could be no n negligible. Wethere- 
fore only measure indices for which iKuntschnerl (120041) has 
provided corrections for these effects (17 indices) though in 
practice we only make use of the R8a, H7A, H/3, Fe5015, 
Mgfe Fe5270 and Fe5335 indices in the following analysis. 
Errors in the line indices were calculated by 500 Monte - 
Carlo simulations of the input spectra. 

Because of the limited number of stars in the Lick li- 
brary observed during the course of the observations, an 
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Table 2. Offsets to the Lick/IDS system derived from Jones 
(1997) library 



Table 3. Major Axis Kinematics 



Index 


Offset 


Number of stars 


trx 
rid A 


— U.oo ± U.Uo A 


lUo 


H<5 F 


-0.16 ±0.03 A 


110 


CNi 


0.003 ±0.002 mag 


117 


CN 2 


0.006 ± 0.002 mag 


115 


Ca4227 


0.00 ±0.02 A 


126 


G4300 


-0.28 ±0.04 A 


128 


H 7A 


0.38 ±0.04 A 


126 


H7f 


0.09 ± 0.02 A 


126 


Fe4383 


0.28 ± 0.05 A 


126 


H/3 


-0.12 ±0.02 A 


128 


Fe5015 


0.23 ± 0.04 A 


126 


Mgft 


-0.08 ±0.02 A 


128 


Fe5270 


-0.07 ±0.02 A 


128 


Fe5335 


-0.04 ±0.03 A 


128 


Fe5406 


-0.06 ±0.02 A 


126 



Notes: Column (1) gives the index name, while column (2) gives 
the mean offset (Lick - Jones) to the Lick/IDS system evaluated 
from the Jones (1997) stars in common with Lick. Column (3) 
shows the number of stars used in the comparison. 



accur ate determination of the offsets (see e.g.. iKuntschnerl 
2000) between line indices measured here and the Lick sys- 
tem could not be determined empiricall y. However, w e made 
use of the stellar library observed by IJonesI <ll997l) to de- 
termine the offsets for indices measured on flux calibrated 
spectra to the original Lick/IDS system assuming that the 
Jones library is well flux calibrated. The spectra were first 
broadened to the Lick/IDS resolution with a wavelength 
dependent Gaussian assuming a spectral resolution of the 
Jones library of 1.8 A (FWHM). Then we compared the in- 
dex measurements in common with the Lick stellar library 
(Worthey et al. 1994) for up to 128 stars depending on the 
availability of the measurements. The offsets and associated 
errors (see Table were derived with a biweight estimator. 
The offsets are generally small (< 0.1 A) but individual in- 
dices can show larger offsets (e.g., H7A, Fe5015). Here, we 
have only considered a single offset per index and ignored 
possible trends with line-strength. Several indices (e.g., H<5a, 
CN2, G4300) show weak evidence for such trends which are, 
however, difficult to quantify (see also Vazdekis 1999). Our 
determination of Lick offsets derived from Jones stars is in 
ex cellent agreement with a n earl ier investigation carried out 
bv IWorthev fc Ottavianl Jl997l Table 9). For the present 
study we apply the offsets listed in Tabic |21 



3 GALAXY KINEMATICS 

Fig 3 shows the result of the kinematical measurements, 
the velocity measurements are based on spectra binned to 
have S/N = 20, the other parameters were derived from 
spectra binned to have S/N = 60. As a check that changes 
in S/N do not affect the measured quantities, the data was 
re-binned to S/N of 30, 40 and 50; the kinematics were then 
remeasured with no significant trends in measured quantities 
being observed. 

We shall discuss each of the parameters measured 
here in turn but in general the measured kinematical 



Radius ["] 


Velocity [km/s] 


<r [km/s] 


H 3 


H 4 


-130.71 


229.33 


106.66 


-0.098 


0.013 


± 


4.15 


5.49 


0.008 


0.007 


-78.66 


245.07 


98.98 


-0.158 


-0.028 


± 


7.9 


8.98 


0.017 


0.014 


-70.52 


244.13 


98.48 


-0.140 


-0.062 


± 


8.04 


10.53 


0.018 


0.017 


-64.41 


242.45 


106.53 


-0.122 


-0.008 


± 


8.29 


9.1 


0.017 


0.013 


-59.47 


242.97 


101.13 


-0.147 


-0.012 


± 


8.68 


8.88 


0.019 


0.013 


-55.11 


242.83 


104.85 


-0.175 


-0.018 


± 


7.97 


9.24 


0.017 


0.014 



Notes: Table continued in electronic format. Tables 4 (minor 
axis kinematics), 5 and 6 (major and minor axis line strengths, 
ages, [ct/Fe] and [Z/H]) are provided in electronic form at 
|http: / /star-www. dur.ac.uk~ dph3man/data. html All data 
provided is binned to S/N = 60. 



data is in very good agree ment with published da t a from 
several authors i ncludi ng: [K^rmenfL^ fc Richst onel ^^j^L 
ICapaccioli et alJ dl993l) . Bender et all Jl994l) and iFisherl 
( 1997) in almost all respects except for the inferred /14 value 
and the rotational velocity at large radii. 

The inner rotation curve of the major axis measured 
here is in good agreemen t between all the data s ets w e 
have examined includin g Illingworth & Schcchtcr (1982), 
Kormen^vfc Richgtona 1^ 921). [Capaccioli^t^flJ 1^1993). 
Bender et alJ Jl994l) . IFisherl dl997l) and lEmsellem et alJ 



ll 19991). At larger rad ii however other authors including 
lOapaccioli et all l)l993l) have measured an essentially flat ro- 
tation curve with rotational velocity around 260fcms _1 . We 
however observe some evidence for some drop off beyond 70 
arcsecs. As observed by several authors, we find no evidence 
for minor axis rotation. Statistics are insufficient at present 
to determine if either of the GC sub-populations can be 
better associated with structures such as discs or spheroids 
within NGC 3115. It is however clear that the bulk of the 
clusters examined rotate in a manner consistent with that 
of the bulk of the stellar content even at larger radii. In fact 
of the 26 clusters with kinematics examined here 21 rotate 
in a prograde manner. 

The velocity dispersion of NGC 3115 is also in good 
agreement with other measured results from the authors 
cited previously. However the data we present here extends 
to significantly larger radii than previous studies. One obvi- 
ous feature of the data presented in Fig. 3 is that the minor 
axis velocity dispersion is considerably higher than that of 
the major axis at equivalent radii. The fact that the mi- 
nor axis data displays a higher velocity dispersion than the 
major axis is not entirely unex pected, as data presented by 
iKormendv fc Richstonel <ll992l) hints at this being the case . 
The 2-D spectroscopy presented in lEmsellem et alJ (Il999t) 
also displays evidence for lines of constant velocity disper- 
sion being elongated in the minor axis direction (at least 
within the inner 5 arcsecs). A difference in measured ve- 
locity dispersion for the major and minor axes can be ex- 
plained if the galaxy consists of a fast rotating, kinematically 
cool disc component and a slower rotating, kinematically hot 
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Figu re 3. Compa rison of measured LOSVD para meters with literature values. Filled circles are from this study, asterisks from 
iBender et alJ ll994l) . and squares from lFisheJ Jl997l) - Triangles show GC data from lKuntschner et alj 120021) de-projected to the major 
axis; filled triangles are the blue sub-population and unfilled the red sub-population. Several GC points are omitted because they exhibit 
counter rotation or have implied rotation velocities of greater than 400kms — 1 . Boot-strapped la error bars are over-plotted for values 
determined in this study, except for major axis velocity and velocity dispersion data where they are omitted in the interests of clarity. 
Typical errors in V and a are < 12 kms" 1 . The major axis data has been mirrored about the minor axis and about the central recessional 
velocity, the minor axis has similarly been mirrored about the major axis. 



spheroidal component. This possibility is re-examined later 
in the light of line index measurements. 



The values of /13 determined here are entirely consistent 
with values determined previously. 



The hi values measured here are generally consistent 
with those measured previously except at large r radii , where 
we observe a value Ka of 0.0, but lBender et all lll994f) find a 
value of around 0.1 this discrepancy could be due to several 
factors including differing experimental set-up and method 
of measuring h3 and h4. 



4 LINE INDICES 



4.1 Abundance Ratios 



In Fig. 4 we s how Mgb vs. <F e> where <Fe> = ((Fe5270 + 
Fe5335) / 2) iGonzaledll993Tl . for the major and minor axes 
of NGC 3115. The GC data points are separated into red 
and blue s ub-populations by (V - I) colour. Model predic- 
tions from lThomas et alJ i2004l) are overplotted for [a/Fe] 
of 0.0, 0.3 and 0.5 with ages = 3 (left most line in each 
group), 5, 8 and 12 (right most line in each group) Gyr and 
have [Z/H] which range from -2.25 (bottom left) to +0.35 
(top right). The effects of age and metallicity are essentially 
degenerate in this diagram, with sensitivity to abundanc e 
ratios maximised. As described in lKuntschner et alJ |_2002), 
abundance ratios are most accurately determined for larger 
ages /metallicities. 
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Mgb [A] Mgb [A] 

Figure 4. Comparison of the [a/Fe] ratios of the galaxy data and GC population through the use of an Mgft vs. <Fe> diagram. In the 
left panel filled and open black circles are binned major and minor axis galax y data from th is paper, filled grey circles show results from 
IFisher et all <1996ft . the open diamond is the central galaxy result from Tra ger et al Fille d black tria ngles are the blue globular 

cluster sample and unfilled triangles are the red globular cluster sample, both frmn lKuTteclTier" et alJ l2002h . Filled black square is the 
error weighted mean of th e blue cluster s ubpopulation and the black square the error weighted mean of the red cluster subpopulation. 
Overplotted are models bv lThomas et alJ feOO.l l2004ll with abundance ratios [a/Fe] = 0.0, 0.3, 0.5, the models have ages 3-12 Gyr and 
mctallicity [Z/H] = -2.25, -1.35, -0.33, 0.0 and +0.35. Right panel shows the unbinned galaxy data presented here, dashed lines denote 
the positions of the binned data. Representative error bars are plotted for spectra located at small, intermediate and large radii, errors 
are a combination of bootstrapped lcr errors and those introduced by a ±5% error in the sky subtraction. 



The measurements for the centre of NGC 3115 from the 
major axis d a ta are in reasonable agreement with those from 
IFisher et all dl996l) and also from lTraeer et all dl99gl) . The 
implied a - element over-abundance of [q/Fe] ~ 0.17 was de- 
termined by chi 2 mi nimization of the data and the models of 
Mgb vs. <Fe> from lThomas et al.1 (120041) . T his chi 2 proce- 
dure, which is similar to the one introduced bv lProctor et alJ 
(2004) fits the set of indices used h ere (H/3, Fe5015, Mg b 
Fe5270 and Fe5335) to the models of iThomas et alJ J2004I) . 
allowing variable ages, [a/Fe] and [Z/H]. A possible break 
in the value of [a/Fe] is visible for Mgb < 4 with an increase 
in [a/Fe] to around 0.25. In comparison, the minor axis data 
starts off at values consistent with the central portion of the 
major axis, and then begins to trend off towards [a/Fe] = 
0.3 much more rapidly. 

A break in the major axis data and differences between 
the two axes can be interpreted as evidence for the existence 
of at least two distinct populations, with typical values of 
[a/Fe] ~ 0.17 and 0.3. Since a move towards lower Mgb and 
<Fe> corresponds to a move to increasing radius it raises 
the possibility of observing radial trends in the strengths of 
other line indices, which will be examined in more detail in 
the next section. 

Our data extends to sufficiently large radii that the 
mean metallicity of the stellar population is similar to that 
of the most metal rich GCs. Fig. 4 shows that on both 
the major and minor axis the mean abundance ratios at the 
largest radii are also consistent with those of the cluster pop- 
ulation although the spread in abundance ratios appears to 



be smaller. No significant population of GCs are found with 
properties similar to the stellar population at intermediate 
radii on the major axis. 

4.2 Radial Profiles of Indices 

Fig. 5. displays the radial profiles of the measured Lick in- 
dices, [a /Fe] and [Z/H]. The valu es of r e used here are those 
listed in ICapaccioli et all il993h for the spheroidal compo- 
nent of the galaxy. 

The differences in [a/Fe] between major and minor axes 
is more clearly demonstrated here, as it is evident that at 
larger radii the major axis data again becomes consistent 
with that of the minor axis. This type of behaviour could 
be understood in terms of changes in the relative contribu- 
tions of disc and spheroidal components, with both major 
and minor axes being dominated at small radii by a nuclear 
component. At intermediate radii the major axis would be 
affected by the influence of the disc component, whereas the 
minor axis would simply be tracing the stellar content of 
the spheroidal component. The change in the behaviour of 
the major axis at large radii could then be understood as 
evidence for truncation of the disc component contribution 
at around 80 arcsecs on the major axis. 

The other metallicity tracing indices (Mgb Fe5270 and 
Fe5335) display similar trends, with slight evidence for 
breaks in the major axis data at around the same radius 
as the one seen i n [q/F e]. This has previously been observed 
bv IFisher et alJ (Il996f) in their edge-on sample of SO galax- 



GMOS Spectroscopy of the SO galaxy NGC 3115 7 



2.2 
2.0 
S 1.8 
1 1.6 
1.4 

1.0 

a o.6 

b. 

O 

ac 

0.0 


-1 

s 

< -2 
«o 
I 

-3 
-0.6 

g -i.o 

-1.8 
-4 

a - 5 



• o o© , 

© 


• © 

© 

© 


4 

• <> 
Jf 


© 


• <5 


; 


*'° : 

: t Jf : 






• »■ — — 

X- - 

© 



3 W 



0.4 

0.3„ 
a; 

r. 

0.2 \ 
o 

0.1" 

0.0 

1.0 



o.o - 

[ 

-0.5 
-1.0 



-2 -1 
log(r/r„) 



-2 -1 
log(r/r„; 



Figure 5. Radial profiles of Lick Indices, [Z/H] and [a/Fe]. 
[Z/H] and [a/Fe] determined by chi 2 minimisation as described 
in Sect. 4.1. Filled black circles show major axis data, open cir- 
cles s how minor ax is data and light grey triangles show data from 
iFisher et alJ il996t) . Error bars are 1-a errors in each bin. Note 
that r e = 93 an d 35 arcsecs fo r the major and minor axes respec- 
tively iCapaccioli et aDll993l) . 



ies (unfortunately they did not examine the minor axis of 
NGC 3115), who found that the Mg2 index is stronger with 
a lower gradient at larger radii on the major axis than on the 
minor axis. As can be seen from Fig. 4 the values for Mgfo 
Fe5270 and Fe5335 determined here for the m ajor axis are 
in goo d agreement with those determined by IFisher et alJ 
(1996). Small differences between the two datasets can be 
attributed to variations in the experimental set-up, the dif- 
ferent methods employed to correct for the broadening effect 
of the LOSVD, and the uncertainty on our correction to the 
Lick system. 



4.3 Age Determinations 

Fig. 6. examines the behaviour of age-sensitive indices along 
the major and minor axes of NGC 3115. Model grids are 
interpolated b e tween the [q /Fe] = 0.0 and 0.3 models from 
iThomas et"al"l (12001 l2004h with [a/Fe] determined as de- 
scribed previously. Note that whilst the [a/Fe] ratios used to 
produce the major axis grids are appropriate for the longslit 
data described here, they are generally not appropriate for 
the GCs which tend to have mean [a/Fe] ~ 0.3. However the 
minor axis grids provide a good approximation to the mean 



[a/Fe] of the GCs and hence the best age determinations 
for the clusters. 

The H/3, H7A and H<5a indices are 
plotted against [MgFe]' where [MgFe]' = 
jM qb x (0.72 x _Fe5270 + 0.28 x Fe5335l iTrhomas et all 
2003). This index was found bv lThomas et alJ (12003^ to be 
independent of [a/Fe] and a good tracer of total metallicity. 

There appears to be a difference in age between the two 
axes with the major axis having a mean age of around 5 - 
8 Gyr and the minor axis an age of around 12 Gyr. This 
would suggest a small amount of star formation may have 
continued in the disc component for several Gyr after the 
formation of the spheroid. 

The age determinations from the H/3 index appear to 
change at lower [Z/H] but in the opposite sense for each axis, 
with the major axis appearing to become younger and the 
minor axis older. As H/3 is relatively unaffected by changes 
in [a/Fe] this cannot be explained as being due to changes in 
[a/Fe] along either axis, but could be explained for the ma- 
jor axis by more recent star formation in the outer parts of 
the disc (in spiral arms perhaps) . Another possibility would 
be that some undetected H/3 emission in the inner regions of 
the disc weakens the observed H/3 index in the inner regions. 
This final possibility however seems unlikely since we find 
negligible signs of [OIII] emission. For the minor axis data 
it would seem plausible that the populations at larger radii 
could be older and would in fact represent the older genera- 
tion of stars also being traced by the GC populations. This 
possibility would seem to gain credence from the fact that 
at larger radii the minor axis data displays [a/Fe], ages and 
[Z/H] values that are entirely consistent with those deter- 
mined from the red GC sub-population. 

The remaining age estimators H7.4 and H<5a also show 
an age offset between the minor and major axis in the same 
sense as H/3, but there is little or no evidence for an age 
gradient. 

Fig. 7. shows the radial profiles of ages determined by 
the chi 2 minimisation described previously. As should be 
expected the trends previously described are obvious with 
the major axis appearing to have an age of 5 - 8 Gyr and 
the minor axis having an age of between 12-14 Gyr. 

A further comment is that despite the use of different 
age sensitive indices our age dete rminations for the GCs 
are entirely consistent with those of lKuntschner et alJ (2002) 
with both GC populations having a mean age of around 12 
Gyr. In fact the agreement is now i mproved as the me a n ages 
of the sub-groups determined bv iKuntschner et alJ (2002) 
varied from ~ 6 Gyr to 12 Gyr, depending on the Balmer 
line being examined. This spread can now be understood as 
being due to the [a/Fe] sensitivity of the H7F and H<5f lines 
used in the Kuntschner analysis, which can no w be corrected 
using the newer SSP models provided by IThomas et alJ 
(2004). 



5 A SIMPLE TWO COMPONENT MODEL 

To test the hypothesis that radial trends in [a/Fe] could be 
explained by intrinsic differences in the disc and spheroid 
components a simple model was constructed. 

The relative contributions of the disc and spheroid com- 
ponents were determined from archived GMOS images of 



8 Mark A. N orris, Ray M. Sharpies and Harald Kuntschner 




-4 - 



0123401234 
[MgFe]' [A] [MgFe]' [A] 

Figure 6 . Age-metallicity diagnostic plots for major and minor axis data. Symbols are as in Figure 4. Overplotted are models by 
iThomas et alJ <200.1 120041) with abundance ratios [a/Fe] indicated in the top right corner of each column, a value of [a/Fe] = 0.25 was 
chosen for the minor axis a s an average of the va lues determined at small and large radii. These [a/Fe] ratios are interpolated from the 
0.0 and 0.3 values listed in lThomas et ail l2004h . the models have metallicity [Z/H] = -2.25 to +0.35, lines indicate ages from top to 
bottom of 3, 5, 8 and 12 Gyr respectively. Error bars for galaxy data are the \-a errors on each of the radial bins. 



NGC 3115 obtained in 2004 (Program ID GS-2004A-Q-9, 
PI R.M. S harpies) in the g , i and r bands. The program 
GALFIT ijPene et, al.ll2002h was used to carry out a simple 
bulge to disc decomposition. Because we only wish to know 
the relative contribution of the two components and we are 
only interested in a simple first order model, a full isophotal 
decomposition was unnecessary and a simpl e de Vaucouleurs 
mode l of the spheroid with r e = 93 arcsecs (ICapaccioli et all 
1993) was subtracted from the original image to isolate the 
disc light contribution. As a test of this approach, we also 
allowed GALFIT to attempt to fit the images with a Ser- 
sic function with variable r e and n, both with and without 
masking of the disc region and with varying starting param- 
eters. The results in terms of r e and n varied considerably 
but the distribution of flux in the residual image remained 
fairly constant, with the disc tending to provide a peak of 
around 30 - 40% of the flux on the major axis. 

The results of this bul ge subtraction can be se en in Fig. 
8. As previously noted bv ICapaccioli et al.l (119881) the disc 
shows considerable flaring in the outer regions, which these 
authors attributed to the disc ceasing to be self gravitating 
in this region. There is also evidence for structure within 
the disc (spiral arms?), which could be taken as evidence 



for some residual star formation events. Note that the inner 
region is not well fit by this model. A more realistic model 
would require several components, but is beyond what is 
required for the present analysis. 

The model can be used to predict [a/Fe] at any point, 
if it is assumed that the residual light traces an enriched 
disc and nuclear component with [a/Fe] ~ 0.0 and that the 
fitted spheroid traces a lower Fe - enriched population of 
[a/Fe] ~ 0.3. By weighting the [a/Fe] value by the relative 
fractions of the two components it is possible to estimate 
the observed [a/Fe] at any point. Fig. 9 shows the result 
of this procedure for both axes. This simple model repro- 
duces the general trends observed with [a/Fe] quickly rising 
on the minor axis and a much more Fe-enriched major axis 
which trends back towards the asymptotic values of the mi- 
nor axis at large radii. The differences between the model 
and observed values can be explained by a number of factors 
including incorrect values for the intrinsic [a/Fe] of the two 
components and the effect of other unaccounted for com- 
ponents. Other factors such as intrinsic gradients in [a/Fe] 
could also play a part. However to first order we believe that 
differences in [a/Fe] between the two axes of NGC 3115 can 
be explained as being due to the existence of at least two 
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Figure 7. Radial profiles of measured age for major and minor 
axes. Ages determined by chi 2 minimisation as described in Sect. 
4.1. Symbols as previously defined. 



distinct stellar populations within the galaxy, with different 
spatial distributions and enrichment histories. 



6 DISCUSSION 

In general all of our kinematical measurements are in good 
agreement with previous studies, the most interesting find- 
ing being that the disc of the galaxy is particularly cold 
(vdisc I &disc > 2) and hence rotationally supported. An 
interesting extension to this work would be to improve the 
statistics of the kinematics beyond 100 arcsecs, which would 
probe the region where substantial flaring of the disc has 
been observed. The fact that the rotation velocity on the 
major axis remains high even at large radii where the inte- 
grated light is dominated by the underlying spheroid implies 
that in NGC 3115 the halo is also rotating significantly (v / 
a ~ 1.5) in the same sense as the disc. 

We have confirmed that the GC population shows clear 
evidence for rotation, in the same sense as that of the disc. 
Statistics on the GC kinematics at present are not sufficient 
to associate particular cluster populations to any specific 
structural feature. 

Our investigation of the Lick absorption line indices has 
shown that the stellar disc component is considerably differ- 
ent from the spheroidal component in almost all measured 
parameters. Most informative in terms of constraints on for- 
mation theories is the implication that the disc of the NGC 
3115 is several Gyr younger, and significantly enriched in 
metals, compared to the spheroid of the galaxy. This has 
been hinted at previously b y studies of the op tical colours 
of NGC 3115, for example ISilva et al.l dl989l) . who noted 
that in B - i the disc was 0.5 mag bluer than the spheroid. 
This colour difference however could be attributed either to 
a younger stellar population in the disc or to the disc hav- 
ing a lower metal abundance compared to the spheroid. We 
have convincingly demonstrated that a lower metal abun- 
dance cannot explain this colour difference since the major 
axis displays a higher metal abundance than that of the 
minor axis. The difference in colour can therefore clearly 
be attributed to an age difference of 5 - 7 Gyr for the two 
components and is also consistent with the observation that 
[cv/Fe] ratios are lower on the inner region of the major axis. 

The red GC sub-population is most consistent with the 
larger radii minor axis data in its measured parameters. This 
suggests that both the spheroid and the red GC popula- 
tion may have formed from the same material at around the 
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Figure 8. Top image shows original GMOS i band image. 
The 2nd image shows the simple de Vaucouleurs model for the 
spheroidal component. The 3rd image shows the residual image 
after the 2nd is subtracted from the 1st. These 3 images are all 
scaled to the peak intensity of the original image for clarity. Con- 
tours show the 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9 fractions. 
The lower image shows the fraction of the total emission provided 
by the residual components. Considerable flaring of the disc is 
visible, as is substructure consistent with spiral arm structures. 
Contours show the 0.1, 0.2, 0.3, 0.4 and 0.5 fractions. 



same time. An interesting extension would be to probe to 
larger radii to see if the observed trends in metallicity and 
[a/Fe] shown in Fig. 5 continue, and eventually lead to the 
minor axis spheroid data tracing stellar populations similar 
to those that make up the lower metallicity blue GC pop- 
ulation. This blue population of GC's could then possibly 
be attributed to an initial burst of star formation during 
the halo formation of the galaxy. One unanswered question 
is whether or not there exists a GC population associated 
with the substantial younger disc component of this galaxy. 
At present the number of GC's for which spectroscopy ex- 
ists is too small to exclude this possibility. Though at present 
none of the GC's examined have line-strength indices com- 
parable to those measured for the disc component. If no such 
population exists, this in itself would prove interesting as it 
would imply that GC's are not necessarily formed whenever 
a sizeable amount of star formation occurs. 
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Figure 9. Model predictions for variation of [a/Fe] for major 
and minor axes. Solid line shows the model prediction for the 
major axis, dashed line displays the model prediction for the mi- 
nor axis, both made assuming [a/Fe] = 0.0 for the disc and 0.3 
for the spheroid. The circles show radially binned [a/Fe] measure- 
ments shown previously in Fig. 5 for the respective axes. Slight 
differences between both sides of the major axis are visible in the 
model predictions, which have been folded about the centre. 

7 CONCLUSIONS 

We have presented new kinematical data and absorption line 
strength measurements for the major and minor axes of the 
SO galaxy NGC 3115, and have compared these measure- 
ments to similar data for the GC system of NGC 3115. Our 
main conclusions are: 

• NGC 3115 has a significant stellar disc component, 
which is both kinematically and chemically distinct from 
the surrounding spheroidal component. 

• The spheroidal component of NGC 3115 is consistent 
with having a uniformly old ~ 10 - 12 Gyr age and [a/Fe] of 
0.2 - 0.3. At large radii the minor axis (which should trace 
the spheroidal component exclusively beyond the central 
few arcsec) is consistent in age, [a/Fe] and metallicity with 
the red GC sub-population, hinting at a common origin for 
the two. 

• The major axis data displays clear evidence for con- 
tamination by a younger (5-8 Gyr old) more chemically 
enriched stellar disc. The observation that the disc of 
NGC 3115 is bluer than the spheroid is primarily an age 
difference not a metallicity difference effect. 

• Previously observed discrepancies in age determination 
between the H/3 and higher order Balmer lines for the GC 
sample can largely be explained by changes in the higher 
order Balmer lines due to varying [a/Fe]. 

• The GC system displays clear evidence for prograde 
rotation in the same sense as the disc and spheroidal com- 
ponents. 
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